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Quantum information processed in strongly correlated states of matter can provide built in hardware protec- 
tion against errors. We may encode information in iiiglily non local degrees of freedom, such as using three 
dimensional spin lattices for subsystem codes or two dimensional spin lattices for topologically ordered surface 
codes and measurement based codes. Recently, in [L. Jiang et al.. Nature Physics 4, 482 (2008)] the authors 
showed how to manipulate these global degrees of freedom using optical lattices coupled to a bosonic degree 
of freedom via a cavity. We elaborate on these ideas and recapitulate two approaches to implement many body 
gates necessary for quantum information processing, both relying on controlled interactions of an ancillary cav- 
ity mode with the spin system and single ancilla particles. The main focus of the present paper is to analyze the 
effect of imperfections such a cavity decay and collective and individual spin decoherence. We present strategies 
to fight decoherence by monitoring cavity decay and show that high gate fidelities can be achieved in the strong 
coupling regime of cavity-QED with state of the art parameters. 

PACS numbers: 03.67.Lx,75.10.Jm,37.10.Jk,37.30.+i 



I. INTRODUCTION 

One way to store quantum information is to prepare it in 
the degenerate ground subspace of a many body Hamiltonian. 
The preparation will be robust if the logical states are stable 
under environmental noise. This can be established if there is 
a gap dE between the code space and its orthocomplement and 
if the logical operators themselves are sufficiently non local in 
nature. In such a case, at sufficiently low temperature T, there 
will be a suppression of individual errors by a factor g^^^/^sr 
and an accumulation of errors that leads to an unrecoverable 
logical eiTor will be unlikely. Models for quantum memo- 
ries and computations based on such codes have been put for- 
ward in [ , , ]. Means to realize the necessary spin lattice 
Hamiltonians, caiTying the protected degenerate ground sub- 
space, exist in systems of interacting atoms or polar molecules 
trapped in optical lattices [ , , ]. Yet manipulating such 
codes is non trivial for the very reason that the states are only 
coupled by global operations. That is, gates on single logical 
qubits are realized by multi-qubit gates, which are notoriously 
hard to implement physically. Other strategies for a digitized 
simulation of many body Hamiltonians use quantum circuits 
or entanglement assistance for teleporation of gates (see e.g. 
Ref. [ ]). However, the quantum circuit approach often takes 
the system outside of the code space mid-circuit and the en- 
tanglement assisted gates require m partite GHZ states to sim- 
ulate a single summand of an m body Hamiltonian. 

In Ref. [ ] the authors addressed this challenge suggesting 
a method to implement the required multi-qubit gates with ul- 
tracold atoms or molecules in an optical lattice embedded in a 
high-quality optical resonator (see Fig. 1). The key idea is to 
achieve controlled, collective interactions of selected subsets 
of spins with the cavity mode interspersed with interactions of 
a single ancilla particle coupled to the same cavity mode along 
with measurements of the ancilla qubit. Using these tools two 
types of gates can be constructed: one is based on the idea 
to teleport a gate from the ancilla qubit to the encoded qubit 



and involves single photon excitation of the cavity mode. The 
other is a geometric phase gate and requires coherent excita- 
tions of the cavity mode along with conditional phase rotation 
controlled by the ancilla qubit. 

While all interactions of atomic (or molecular) spins and 
the cavity field are obtained in the dispersive limit, cavity de- 
cay and spontaneous emission will necessaiily affect the gate 
fidelities. In the present paper we complement our recent pro- 
posal [ ] and provide a detailed analysis of the fidelity of gate 
operations in the presence of decoherence. We identify the 
requirements on physical parameters of the system and show 
that high gate fidelities can be achieved in the strong coupling 
regime of cavity-QED with state of the art parameters. Sig- 
nificant experimental progress has been made toward coherent 
control in this regime [ , , ]. We also identify some adap- 
tive protocols, based on counting of photons leaking out of the 
cavity, which can improve the gate fidelity in the presence of 
cavity decay. 

The paper is organized as follows. In Section II (supple- 
mented by Appendix ) we introduce examples of protected 
quantum memories and describe their possible implementa- 
tions in optical lattices. In Section III we present the two 
gate operations, the single photon protocol and the geometric 
phase gate. Section IV (supplemented by Appendix ) provides 
a comprehensive treatment of the impact of cavity decay on 
both protocols along with a discussion of strategies based on 
monitoring the cavity decay. Section V finally treats other de- 
coherence effects, in particular individual and collective spin 
decay. We conclude with a summary of the results. 



II. PROTECTED QUANTUM MEMORIES 

We begin with a brief review of the basic properties of 
ground state quantum memories in three models: subsystem 
codes, surface codes, and spin— 1 chains used for measure- 
ment based computation. 



2 




FIG. 1: Schematic depicting a step in implementing a many body 
operation on trapped atoms or molecular spins mediated by a cavity. 
The lasers for optical trapping are indicated as is the region of max- 
imum cavity field strength. Planes of spins can be moved into the 
strong cavity coupling region by by stretching the lattice along one 
dimension (here along z) by e.g. rotating the z trapping beams away 
from the z axis by an angle ^. The lattice beams are assumed tuned 
far from cavity resonance so that the mirrors are transparent at that 
frequency. One can also use additional lasers with shaped intensity 
profiles [2 1 , 22] to make a particular region C of the lattice optically 
active such that only spins in that region can interact with the cavity 
mode. 

A. Subsystem codes 

We focus on a specific subsystem code, the three dimen- 
sional spin lattice model studied by Bacon [3] also known as 
the 3D compass model. The system is comprised of qubits 
residing on the vertices of an n x n x « simple cubic lattice (n 
odd), as illustrated in Fig. 2. The interaction Hamiltonian is: 

^cp = -/( L ctV + L a>o> + ^ <fa'), (1) 

X— links y— links z— links 

where J > Q. We label the single qubit operators afj f. for 
o" at site x — i,y = j,z = k. It is unclear how the gap of 
this system scales with size though numerical evidence for the 
2D compass model suggests that the gap scales as 2^" [12]. 
Nevertheless, one could encode in a small system to get the 
benefit of ground state protection. 

Stabilizer operators are generated by 2 (n — 1) operators as- 
sociated with adjacent planes, including operators in the 
x — y plane and o' operators in the y — z plane: 

^ = n <j,K^+i.j,k^ ^ = n '^,jj<^i,.k+v 
j.k=\ i.j=i 

We can encode one qubit of information in the ground states 
of Hep, with logical operators X = YTj k=i'^l j k (product of 
operators in the y — z plane) and Z = 11?;= i ^; ; i' product of 
operators in the x — y plane (see Fig. 2). Properties of this 
subsystem code are discussed in Appendix A. 



B. Surface codes 

Consider the 2D lattice where each edge of the lattice repre- 
sents the location of a spin- 1/2 particle with a coupling graph 
such that particles on edges which meet a vertex interact via 
Hy = nygstarjvjCTy and edges which surround a face interact via 
Hf = Yljfz^fO'j. The operators and Hf always collide on 
an even number of edges and hence commute. Fuithermore, 
they assume eigenvalues ±1 and information can be encoded 
in the +1 coeigenspace of these operators. By choosing the 
so-Ccilled surface-code Hamiltonian [ , ]: 

H,^rf = -UY,H,-jJ^Hf, (2) 

f 

with U,J > the code space corresponds to the ground 
subspace with degeneracy that depends on the topology: 
dim!J-(gr = 2^*+'' where g is the genus of the surface and li is 
the number of holes [ ]. Designing lattices with genus g >0, 
e.g., on the surface of a torus, would be challenging, but al- 
ternatively it would be possible to create several holes (h > 0) 
in a planar lattice with boundary by for instance deactivating 
regions of the lattice with focused far detuned lasers. Further- 
more, the planar code with rough boundaries can also provide 
a twofold ground state degeneracy of [ ]. Since they are in- 
sensitive to local perturbations, the degenerate ground states 
provide a good quantum memory. A caveat is that thermal 
fragility of the topological order implies a threshold temper- 
ature Tihres above which logical information decoheres. This 
temperature scales as Tlhresh ^ Sg/lnn [15] where n is the lin- 
ear dimension. 

The code states are coupled by the logical operators: Z = 
Ylj^f^^Op and X = Ylj^f^j^a-^j where the configurations Cz{Cx) 
are strings on the lattice (dual lattice) as illustrated in Fig. 3. A 
more realistic two body Hamiltonian described by anisotropic 
nearest neighbor Ising like interactions on honeycomb cou- 
pling graph was proposed by Kitaev [ ] which in a certain 
parameter regime yields an effective Hamiltonian unitarily 
equivalent to //surf. 

C. Ground code computing 

Another means of processing information in many body 
states is via ground code measurement based quantum compu- 
tation (GMQC). Here the information is processed in degener- 
ate ground states of a gapped Hamiltonian over a two dimen- 
sional lattice of spins, and computation flows by sequential 
measurement on the constituent spins. It was shown in Refs. 
[16, 17] that it suffices to have a nearest neighbor only inter- 
action between spin— 1 particles to realize GMQC with single 
and two spin measurements. In the protocol of [17], each log- 
ical qubit is stored and processsing in the ground states of spin 
chain with spin— 1 particles in the bulk and spin— 1 /2 particles 
on the boundaries interacting via the so-called AKLT Hamil- 
tonian [IS] 

Haklt =J[Y, +^0,1 +Pn,n+i] 
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with 7 > 0. Here Pj j^i is the projector onto the spin-5 ir- 
reducible representation of the total spin for particles j and 

1, i.e. = i (S, • S,+i + i (S; • S;+i )2) + ^ and P^/? = 

|(l6 + Sj -Sji), where S,s are spin-1, 1/2 representations of 
5m(2). The ground state of //aklt is non-degenerate but after 

3 /2 

turning off the interaction Pq\ and measuring the qubit at lo- 
cation 0, the chain is initialized into a logical state of a two 

fold degenerate ground subspace of H'^i^^j = Ha^lt — Pq^i- 
Computation flows by measuring spins along the the chain 
and pairs of neighboring spins in parallel chains with the final 
output state on the last qubit located at + 1 . Such a model 
is protected by the hardware from errors but for fault toler- 
ant computation it may be useful to employ multiple chains 
for each logical qubit using a quantum error correction pro- 
tocol where the state of the last qubit of one chain is tele- 
ported into the ground states of a freshly prepared chain. The 
ground states of H'^j^j^j are connected by the string operators 
2:A' = e"^^*=i'^*®0^^j. Hence, given the string initialized in 
the state induced by measuring the qubit at position 0^ in |J,Bq) 
we can teleport the logical state of qubit A into B via the uni- 
tary f/ = |0) (0| (g) 1 + 1 1) (1 1 (g)E-"(;)- We describe in Sec. UI C 
how this could be done by mapping the state of the last qubit 
of one chain to a photon, then allowing the photon to interact 
with the new chain to generate the many body operator U. 




FIG. 2: A cubic lattice with anisotropic nearest neighbor interaction 
encoding one logical qubit in physical qubits. Spins interact ac- 
cording to a Hamiltonian Hep with o"ct" interactions among each 
pair of spins in directions a = x.y.z- This sytem can be implemented 
e.g. using microwave induced dipole-dipole interactions between po- 
lar molecules as shown in [ ]. The planar operators and Lf which 
involve a product of and CT~ operations along the respective y~z 
and x — y planes act as logical X and Z operations on the code. 



E. Physical implementation 



D. Universal operations 

In the remainder of this paper we consider how to generate 
unitary evolution by the many body operators 

SI=U4 (3) 

where C is some set of spins and \C\ = m. We may be in- 
terested in performing string operators on surface codes [ ] in 
which case C is a set which defines a connected string of spins 
in a 2D lattice (see Fig. 3), or perhaps we wish wish to per- 
form encoded operations on the [[ir' , 1,«]] subsystem code in 
which case is a plane of spins (see Fig. 2). One method uses 
a single photon to generate the many body interaction and the 
other uses a geometric phase gate. 

Given the abiUty to generate arbitrary rotations e^^^^ , e'^^^c 
and the CNOT gate, and measurements of Z, universal quan- 
tum computation can be achieved. For both surface codes 
and subsystem codes discussed above, the logical CNOT op- 
eration can be done transversally between two code blocks. 
The logical operator Z (or X) is a global operator that cou- 
ples multiple spins over the encoding block. The measure- 
ment of Z can be achieved by measuring individual physical 
spins. However, it is a non-trivial task to perform the logical 
rotations e"l''^r. The center aim of this paper is to provide two 
approaches to perform such logical rotations. 



Methods for analogue simulation of the three models above 
have been proposed using trapped atoms [4, 19, 20] or 
molecules [ , ] in optical lattices. For trapping of polar 
molecules with lattice spacings of ^ 200 nm, the interaction 
strength using microwave induced dipole-dipole interactions 
can be as high as 7 ~ 27t x 10 kHz with decoherence dom- 
inated by spontaneous emission of optical photons at a rate 
Y ~ 27t X 1 Hz. Such small lattice spacings are possible us- 
ing trapping lasers tuned to the properly chosen molecular 
transitions [ ]. The implementation //imp will not be exact. 
There will be spurious longer ranger interactions and devia- 
tions from the required symmetry of the nearest neighbor in- 
teraction. Yet such deviations need not break the code. For 
example, the latter two models //surf and //aklt have a gap 
which provides for some resilience to imperfect implementa- 
tion. Furthermore, for an implementation of //^p with polar 
molecules, the microwave fields that induce the interactions 
are linearly polarized and the errors are the form of products 
of pairs of Pauli operators which hence preserve the time re- 
versal symmetry of the model. For this reason, the degeneracy 
of the code states is preserved and the gap condition for suffi- 
ciently small systems, can be maintained provided the devia- 
tions are small. Just how small depends on a detailed compu- 
tation of the energy gap to excited states. By optimizing the 
microwave beams that induce the dipole-dipole interaction, its 
found that for a 4 spin configuration on a 3D trine, the devi- 
ation of the implemented Hamiltonian //imp to the target //cp 
is: ||//imp — //cpll < 10^'*7, where the norm is defined as the 
supremum norm of the traceless part of the operator [: ]. 
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FIG. 3: A planar code which encodes one logical qubit in the ground 
states. There is a spin-1/2 particle (filled dot) for each edge of the lat- 
tice. The interactions of the local Hamiltonian //surf ars along edges 
that bound a face /, and edges that meet at a vertex v. The strings 
Cx,z indicate paths of products of operators that are logical op- 
erators on the code. 

III. IMPLEMENTATIONS 

We focus on physical systems in which an optical lattice is 
placed within a high-finesse optical cavity. We assume that 
we can control the coupling between the cavity mode and a 
selected set of spins in the optical lattice. Such selective ma- 
nipulation can be achieved using a control laser beam with ap- 
propriately shaped intensity profile [8, 21, 22]. Alternatively, 
if selected spins form a simple pattern, such as occupying all 
the sites in a straight line or a plane, we can stretch the lattice 
[23] and couple the cavity mode with the spins only located at 
certain line or plane. 

We assume a spin dependent dispersive coupling of our lat- 
tice to a cavity mode: 

Vz^gzfl^fl^o}, (4) 

where gz is the dispersive coupling strength, and C is the set 
of selected spins that can be associated with string operator for 
the toric code or AKLT model, or planar operator for the sub- 
system code. In addition, it will be convenient to have a spin 
dependent coupling of a single ancillary spin to the cavity: 

VA^gAa^a\\)A{l\A. (5) 

with coupling strength gA, which can be activated by bringing 
an ancillary particle into contact with the cavity mode, allow- 
ing them to interact from some time (without interacting with 
the spin degrees of freedom of the system particles), and then 
de-activated by removing the ancillary particle. 

The dispersive coupling of the atoms or molecules to the 
cavity can be achieved via a state dependent AC stark shift 
[Fig. 4(a)]. When the cavity frequency (0^ is detuned by A far 
off resonant from the excited states, then there will be a dis- 
persive interaction that is photon number conserving. These 



excited states could be electronic excited states in the case 
of an optical cavity or rotational excited states in the case of 
a microwave cavity. If, for example, the photonic mode is 
0+ polarized then there will be a differential shift on the |0) 
and |1) spin states of the polar molecules due to the different 
angular momentum coupling coefficients co.i for the ground 
to excited state transitions [ ]. Up to a constant this inter- 
action is equivalent to Vz where gz ~ (cg — c\) 2va'' ^^^^ 
V the effective mode volume of the cavity and d the optical 
dipole moment of the polar molecule. The ancillary particle, 
with a different state space such that only one state 1 1)^ inter- 
acted with the cavity mode could be brought into and out of 
the cavity with optical tweezers to generate evolution by Va- 
Recently, a dispersive interaction between atoms in an optical 
lattice and an optical cavity was proposed as a way to measure 
quantum phases [ ]. There the cavity frequency is chosen far 
off atomic resonance such that the interaction is of the type in 
Vz except that there is no spin dependence of the atoms so 
that the atomic number operator rather than is measured. 
When coupling via an optical cavity it wiU be important to 
choose the lattice spacing along the cavity axis commensurate 
with the cavity spatial mode spacing in order to ensure equal 
coupling to all spins. 

It would be advantageous to have a way to turn on and off 
the coupling between the spins and the cavity. This is pos- 
sible using cavity assisted Raman pulses. Here the idea is to 
introduce an auxiliary classical field £2 detuned by A from a 
transition g ^ e and have the cavity field detuned by A + 5 
from a different transition gl e. For A ^ £2,y,g where y 
is the linewidth of the excited state, then there is an AC stark 
with ^eff = ii^(l/A+ l/(A+5))/2, and the auxil- 
iary field can turn the coupling on and off. It is difficult to 
find such a closed optical transition in polar molecules ow- 
ing to the complex state space which tends to couple ladders 
of vibrational levels; although nearly closed transitions do ex- 
ist (see [ ]). It is possible to find microwave cavity assisted 
Raman processes for example by choosing the cavity field to 
be tuned near the N — 1 N = 2 transition and the auxiliary 
classical field on the N = ^ N ~2. The latter would have to 
be a strong field to couple via quadrapole transitions. Depend- 
ing on the intermolecular spacing, it may be necessary to use 
at least two cavity frequences to obtain only single spin inter- 
actions without inducing spurious dipole-dipole interactions 
(see methods of [ ]). 



A. Single photon protocol 

Many body gates can be generated with a single cavity ex- 
citation coupling to the lattice. The basic idea is to teleport the 
quantum gate from the probe qubit to the encoded qubit of the 
lattice spins. One choice for the probe qubit can be the photon 
number states of the cavity mode, with zero or one excitation. 
Alternatively, we may introduce an ancilla spin as our probe 
qubit, which couples to the lattice spins via the common cav- 
ity mode. In the rest of the discussion, we will use the second 
choice, with the advantage that it is more convenient to ma- 
nipulate the ancilla spin compared to the photon number states 
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FIG. 4: Cavity-assisted gate based on single photon approach, (a) 
The energy levels of a selected memory spin (|0) and 1 1)) interacting 
dispersively with the cavity mode, which implements the Hamilto- 
nian with spin dependent dispersive interaction. The coupling coeffi- 
cient is gz = g^/\ with single-photon Rabi frequency g and detun- 
ing A from the excited state \e). (b) The energy levels of the ancilla 
spin (different from memory spins) and the cavity mode for the sin- 
gle photon approach. A different control laser with Rabi frequency 
O.A (t) connects the states 1 1)^^ ® |vac) and 1 1')^ ® a' |vac), and en- 
ables coherent creation and absorption of a cavity photon conditioned 
on the ancilla spin, (c) Cartoon illustration of the procedure for the 
implementation of single-photon approach for controlled many body 
gate: (1) Initialize the ancilla spin (the left highlighted spin) in a su- 
perposition state ot |0)^ + P 1 1)^ (blue for |0)^ and red for 1 1)^), with 
no photon in the cavity and state |\|/)^ for the topological memory. 
(2) Coherently create a cavity photon (orange shade) for ancilla spin 
state |I)^ (upper branch); no photon is created for ancilla spin state 
|0)^ (lower branch). (3) Switch on the interaction between the cav- 
ity photon and the selected spins. If there is a cavity photon (orange 
shade), a non-trivial evolution (pink dots) is implemented. (4) 
Turn off the interaction and coherently absorb the cavity photon into 
the ancilla spin. Finally the state a |0)^ ® Ixif)^ + (3 1 1)^ ® 5^ is 
prepared. 



of the cavity mode. 

The procedure of teleporting the quantum gate is summa- 
rized in Fig. 5. First, we couple the probe qubit and the lattice 
spins. Then, we perform the quantum gate over the probe 
qubit. Finally, we measure the probe qubit, which determines 
the Pauli frame of the encoded qubit of the lattice spins. Af- 
ter these operations, the many body gate has been effectively 
applied to the lattice spins. 

Before we give the procedure of coupling the ancilla 
spin and the lattice spins, we first describe a technique 
that allows use to entangle the ancilla spin and the cav- 
ity mode, using the energy levels shown in Fig. 4(b). 
Suppose the ancilla spin and cavity mode starts with 
state |1)^ (8) I vac). When we slowly increase the Rabi 
frequency of the control laser i2(f), the ancilla spin and 
cavity mode adiabatically follow the dark state |A(f)) = 



1+) — 



FIG. 5: Gate teleportation circuit for arbitrary x-rotation Xq = e'^^c 
on the memory. The circuit represents the following procedure: 
(1) use the CNOT gate A[x] to entangle the probe qubit (upper 
line) and the memory (lower line with a slash), (2) projectively 
measure the probe qubit in a rotated basis, and (3) perform an en- 
coded Pauli X gate over the memory conditioned on the measure- 
ment outcome. An analogous circuit (using A [z] and classically 

controlled Pauli Z gate) implements Zg = e'^^c , hence by Euler de- 
composition, any gate can be formed on the memory using three 
controlled operations. For the geometric phase gate scheme, we 
can actually implement rotations of the encoded qubit without the 
probe qubit, e.g., x-rotation of the encoded qubit can be decom- 
posed as e'^^ =£)(-|3)£) 

|ap| = e. 



ae'i-^) £)(P)D (ae'i^), and choosing 



|vac)-£2(f)|l%®flt|vac))/y^ 



where is the single-photon Rabi frequency. Since |A(f)) 
approaches 1 1')^ (g) |vac) for Q.{t) ^ g', we effectively 
transfer the state |l)^(g)|vac) to — |1')^ (g) a''' |vac). Mean- 
while, nothing happens if the initial state is |0)^ |vac). 
Since adiabatic state transfer is a coherent process, the relative 
phase between 10)^0 |vac) and |A(f)) is maintained. Thus, 
we create an entangled state bewteen the ancilla spin and the 
cavity mode 



1 



V2 



(|0)^(K)|vac) - |l')^(g)fl'''|vac)). 



Similarly, we can reverse the state transfer (from — 11')^ (g) 
fl^ |vac) back to \1a) (8) |vac)) by adiabatically decreasing the 
Rabi frequency i2(f). 

Using the technique described above, we can implement the 
controlled many-body operations [Fig. 4(c)] as the following: 

• The ancilla spin (probe qubit) starts with state a |0)^ + 
P|l)^, the cavity mode has no photon |vac), and the 
topological memory is in state l'^)^- 

• Then we adiabatically turn on the control laser Rabi 
frequency £i(f) and coherently create a cavity photon 
if the ancilla spin starts in |1)^, so we obtain the state 

a |0)^ «) |vac) - p 1 1% (K) a''" |vac) . 

• Next we switch on the interaction between the cavity 
photon and selected spins for time T, during which the 
topological memory undergoes an evolution S'^ if there 
is one photon in the cavity (i.e. |vac)). 

• Finally, we switch off the interaction and adiabatically 
transfer — \l')y^(ga' |vac) back to jl)^® |vac). 

After these four steps, the cavity mode restores the initial 
state I vac), while the ancilla spin and the topological mem- 
ory evolve from (a |0)^ + p |1)^) «) |\|/)5 to a|0)^ |i(/)5 + 
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P |1)a ® ("O^-^c which is the controlled operation up 
to a known phase. 

By choosing S'^^ as the encoded Z operator, we implement 
the controUed-Z gate between the ancilla and the lattice spins, 
which can be converted into the CNOT gate by conjugating all 
lattice spins with the Hadamard gate. Therefore, we can im- 
plement all the quantum gates appeared in the circuit in Fig. 5 
and achieve unitary evolution of many body operators. 



B. Geometric phase gate 

An alternative to using a single Fock excitation of the cavity 
mode is to use coherent state control to perform a geometric 
phase gate [26]. Since these operations use only linear opti- 
cal elements, they may be much easier to realize in experi- 
ment. The mechanism makes use of two basic operators, the 
displacement operator D(a) — g""^-" " and the rotation op- 
erator R{d) = e'^"*" which satisfy the relations: D(p)D(a) = 
e'3(P«*)/2D(a + P), and R{Q)D{a)R{~Q) = D{ae'^). Putting 
these primitives together, one can realize an evolution 

g-iH^t ^ D{-^B)R{QC)D{-a)R{-QC) 

xD{^B)R{dC)D{a)R{-eC), (6) 



• Perform the sequence of steps: 



according to an effective Hamiltonian 

//intf^ |aP|sin(9C + (j)), 



(7) 



where (|) = arg(a) — arg(P) and C is an arbitrary operator com- 
muting with cavity operators a and a ' . Picking in particular 

9 = k/2 and C — Y.jec'^j^ which requires driving of the cav- 
ity with coherent fields and interactions of the cavity with the 
spin system according to Eq. 4, the simulated Hamiltonian is 



Haat = sin((|) + m7i/2) | aP | ]~[ a] 



Hence by choosing C = L/'a-i j k ^l' = 0, for the case 

r 

m odd, we can simulate evolution generated by S^. What is 
required is a spin dependent coupling between each qubit on 
the plane and the bosonic channel. 

A non destructive measurement of the many body opera- 
tor S'^ is possible using the assistance of an ancillary particle 
A that also couples the bosonic channel in a spin dependent 
manner. Controlled displacements can be implemented by the 
sequence 

|O)A(O|0l + |l)A(l|®£>(P) = £»(|3/2)7;(7tfl^fl|l>^(l|) 

xD{-^/2)R{-na^a\l)A{l\). (9) 

The protocol for measurement of (5'^) is as follows 

• Prepare the ancillary particle in the state \+y)A = 
(|0)a + /|1)a) /V2, and the field state in the vacuum. 



(8) 



„-it6VA„-it5Vz 



f/ =D(-p/2)e-'"^'»D(p/2)e-"''VAg 

which returns the cavity to the vacuum state. Choose 
parameters satisfying: -gAt2 = gAti -gA^e = SaIi 
^, ~gzh ^ gzk -gzts = 7i/2, and lap I I mod 2k. 
The displacement operators are generated by Vb and we 
pick the relative phase (|) = arga — argp — 7t/2(0) for 
the number of spins m— \C\ even(odd). Reversed evo- 
lution during the interaction time steps f i , f2 , , f6 can be 
achieved by setting gz^A —gZA^ ^.g. by changing the 
sign of the detuning. 

Assuming the system state was initially \\\t)s the joint 
state of system and ancilla is now 



-y)A = 



\yif)s\0)A+f{m)S'M)s\l) 



V2 



where f{m) = (—1)^ for m even and f{m) = ( — 
for m odd. 

• Measure the ancilla in the basis \±x)a — (|0)a ± 
\1)a)/V2- The probabilities for the measurement out- 
come Sx~±l are 

P{s., = ±l) = \il±f{m){Sl)). 

This protocol can be used to measure the product operators 
5^/5^, where C, C' are possibly overlapping configurations of 
spins, via the iterated sequence: [Yliec' ^iWi^ieC ^'W ^i*-^ 
Hi — (}jg Hadamard gate on qubit /. Measure- 

ment errors can be ameliorated by redundifying the state of 
the ancilla before measurement using many, possible faulty, 
CNOT gates offline between that ancilla and many others pre- 
pared in state |0). 



C. Gates on spin 1 particles 

We conclude this section with a brief discussion of generat- 
ing many body gates for teleportation of quantum information 
in spin— 1 chains as prefaced in Sec. 11 C above. This could 
be done in the context of a lattice embedded in a cavity by 
mapping the state of the qubit located at position + 1 of 
logical chain A to the state of a photon and using the same 
procedure as described above. The goal is to then gener- 
ate the unitary operator U = |0) (0| (g) 1 + 1 1) (1 1 (g) E', where 



0^_i_i, on a new chain B initialized in state 



IIbq)- This is generated by an identical procedure to that for 
controlled string operators on qubits but instead of Eq. 4, we 
use 

Vx = gza'aif^ \S^ = 0),(5-'^ = 0| +0^+i), (10) 

7=1 
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where gz is the dispersive coupHng strength. This kind of state 
dependent interaction between a photon and spin— 1 particles 
can be reahzed using e.g. polarization to differentially couple 
the internal states of the particles. 



IV. GATE FIDELITY WITH CAVITY DECAY 

Cavity field decay at a rate k acts as a source of error for the 
many body interactions which it mediates. For the protocols 
above where the system of spins interact with the cavity field, 
the joint state can be decomposed as 

P(0= E (^AJMjA',J',M'j(\^'J^Mj){A',J',M'j\ 

AJMj,A'J',M'j 



®|OCM,)(pM;lj(f), (11) 

where A,J,Mj are labels corresponding to A labeled irreps 
with total angular momentum J and J' projection Mj (see 
[27]), andaA,y,M,.A',/'.M; =Tr[p(0)|A',/,M;)(A,y,My|]. The 

field states |0Ca/j), IPa/j), may depend on the angular momen- 
tum projections. We describe evolution for the case where the 
field states are Fock states for the single photon probe proto- 
col and when they are coherent states for the geometric phase 
gate. Depending on the protocol used, the field states them- 
selves may be entangled with the state of an ancillary spin but 
we focus on computing fidelities for evolution steps where the 
ancilla is non interacting. 

Consider the evolution during an atom field coupling stage. 
The equation of motion for the joint state is 

= -/[Vz,p(f)] + -i2ap{t)a'^ -a^ap{t) - p{t)a'^a). 

(12) 

The evolution conserves the quantum numbers A, J hence we 
can compute the action e^'A^-'-^j{Q) where: 

A^^^Mj{t)^\A,J,Mj){A',J',M'j\{t)^\aMj){K'J{t)- 
The solutions are easily verified to be given by 

~ b'L »./ (t) 



,2=0 



X fl"A^-'''^^(0)(fl^)"e('2i'z<-V2)«+«' (13) 



where 



^ e-[^+i2gz{Mj-M'j)]t 

^^J^'S^^ = K+i2gz{Mj-M'j) 



(14) 



The evolved state is then 



In order to evaluate the performance of many body oper- 
ations mediated by the cavity we will calculate the process 
fidelity for implementing a many body gate U = e^'^^'c where 
= riyGC^;' ™ ^ configuration C of m — \ C\ spins. Many 
body measurements of the type (5"^) are obtained using the 
many body gates with rotation angle X = %I2 as a primitive. 
The details are given in the following subsections with the 
main result that the expected error scales like k/|§z|- 



A. Single photon mediated gate 

Consider a protocol where we begin with the separable 
state; 

p(0) =p5(0)®pf(0) 

KJMj,/\' ,J' M'j 



where 



A^-"<(0) = |A,/,M,)(A',/,M;|®be,=|)(ye,=|l, 



and the field states are defined 



bej =cose+|0)+sine+|l), be 



be 



Here the states |0) , 1 1) could denote photon number in a given 
mode or a single photon in two orthonormal modes where 
only mode 1 1) interacts with the spins. 

The discussion here can be generalized to the case that the 
cavity mode is also entangled with an ancilla spin, with the 
mapping from { |0) , 1 1) } to { |0)^ ® |vac) , - 1 1')^ (g) |vac) } . 
In principle, the mapping does not hold for the cavity decay, 
because the decay process without the ancilla spin (from |1) 
to |0)) cannot be mapped to the decay process with an entan- 
gled ancilla (from 1 1')^ ® |vac) to 1 1')^ ® |vac)). However, 
if we are only interested in the fidelity of the lattice spins, the 
analysis here is still valid, because the reduced density matri- 
ces for the lattice spins (after tracing out the cavity mode and 
the ancilla spin) are the same for both cases. 

For atom field coupling over a period t, we have from Eq. 

13 



A^-^^(f) = \A,J,Mj){^,j',M'j\®-{{l + b^^,^,{t)m{0\ 

_^^{i2gzM'j-mt |0) (1 1 + e(-'2«z'W^-V2>| 1) (0| 
^g(-'2^;z(My-M;)-K)f|^^I^J| 



If we choose the interaction time gzl = f , then the action on 
the state is 

p(T) = e^MO) 

= ^(ps(O)®|O)(O|+P(ps(O))0|O)(O| 

+e-'^^/2p5(0)e''^-''®|0)(l| 
+e-W2e-''^^p^(O)0|l)(O| 



+e-'^e-'"-'>s(0)e"^-'^®|l)(l| 



(15) 
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where 



with 



AJ,Mj,A'J'.M'j 

x\A,J,Mj){A',j',M'j\. 

An ideal many body gate resuhs if K = and we measure the 
photon in a rotated basis. Consider the case where the number 
of spins, m is odd. If we measure the photon in the basis |ye^), 
then we obtain outcome ±1 with equal probabilities p± — ^ 
and the resultant state is 



Vmj.m:, = (cosx + sinxe-™"0(cosz + sinxe'''<)). 

Whereas the action of the map £ on the same basis by multi- 
plication by a complex number: 

'E(Ps(0)) = X! ^KJMj.MJ'M'jQmjM'j 
A,J,Mj.A',J'.M'j 

x\A,J,Mj){A',J\M'jI 

where if we fix 9+ = X such that we approximate the unitary 
U then 



, _ Tr[p(T)|9±)(9±| 
''^^^ " Tr[.] 



= e-'9±(-')"'"'^~cp,(0)e'8±(-')"'-'5h 



Say the target evolution operator is e^'^+(^') but we ob- 
tain the measurement result — 1 . Such an outcome is corrected 
for applying the locally generated unitary 5^. The case where 
m is even is handled in a similar way but we measure the pho- 
ton in the basis 

|jc9+) =cos9+|0) + ;sin9+|l), |xe_) = -/|ji:e^+„/2). 

The measurement outcomes ±1 are again equiprobable and 
the conditional state is pf (t) = e''^^'"'^cps{0)e'^^'"'^c. The 
correction procedure given the outcome — 1 is as before. 

We denote E the map for imperfect implementation of the 
gate U — e^'^^c in the case of nonzero K. For our conditional 
protocol this can be represented as: 

'E{psm^p+S+ipsm+P-S'cS-{Ps{0))S'c, (16) 

where 

S±{Ps{0)) = ^ (x{Q±)psm{Q±y+P{ps{Q)) . 

We assume the product of local unitaries 5*^ can be imple- 
mented with perfect fidelity. 

For the remainder of this subsection, we fix m odd and 
m mod 4=1. The other cases follow in a straightforward 
manner. From Eq. 1 5 we have 

X(9±) = cos9±l + sin9±e-''"/V'''-'\ 

Notice that now the probabilities for measuring the photon in 
|9±) are not necessarily equal, rather one finds 

P+-P- = {l-e-ncosi2Q+). 

The process fidelity fpjo(E, f/) measures how close a quantum 
operation E is to the ideal operation U as measured by some 
suitable metric. The fidelity measure we use is the overlap 
between the induced Jamiolkowski-Choi state representations 
of the operations [28]. The action of the unitary on a complete 
operator basis {\A,J,Mj){A',J',Mj\} is multiplication by a 
unimodular number: 

U\A,J,Mj){A',j',M'j\U' =VMj^!^,\A,J,Mj){A',fMjl 



QmjM'j 



1 r 
2 

X 



cos2x(l +%,,m;(t)) + sinxcosxe"'"'/'*'''''^! 

+e-™<''^-''^Hsin'x(l+V.M;W)-sinxcosx 
X (e-''^^^ + e™'^^ ) + cos^ x 



Of course for the case of no decay, K = 0, then Qf^j j^^t = 

The process fidelity is readily computed using the fact that 
the noise map 'E(pj(0)) commutes with the target unitary U . 
Hence, we can compute the fidelity which measures how close 
the noisy map 'E'(ps(0)) = U^'E{ps{Q))U is to the ideal op- 
eration, i.e. the identity operation: 

Fp,o(£,f/) =Fp,o(£',/) ^s.s' («>+|p£'|1'+)s,s'- 

Here we are computing the overlap of the Jamiolkowski-Choi 
representations of the maps as states in the Hilbert space 
9{s ® ^s' containing our system space and a copy each with 
dimension D: 



1 



^ \A,JMj)s®\KJ,Mj)s'. 



AJMj 

p^, = Is®'E's,{\^+)s,s') 
_ 1 
" ~D 



MjM'j^MjM'j 



AJ Mj ,A' J' M'j 

X |A,7,My)s(A',/,M;| ® \A,J,Mj)s'{A',j',M[ 



Hence 



A.A! J.J'Mj=-Jm'j=-J' 



where 



MjM'j 



Now RmjMj — 0> and the for off diagonal elements, for 
'^l\gz\ < 1 we find. 
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2 2\gz\ ^ ^yigz' 



1- 



(_l)(2My + l)/2 3i„(4^) ' 
2 2{Mj-M'j) j 2\gz\ 



(-l)PMy + l)/2sin(4)c) , cos(4x) + l 
A{Mj-M'j) ^ 2(Mj-M',)2 



r 



Mj—M'j even 

.^(3+cos(4x)) -)(_^)2 M,-M;odd. 

02) 



The above expression is made a bit simpler by counting c'J, 
the number of inequivalent spin— 7 irreps of the m fold sym- 
metrized direct product of SU{2), where m is the number of 
spin— 1/2 particles in the system. The dimension can be com- 
puting using Young tableau [ ]: 



2J+1 



m 

- J 



Making use of the fact that Qm'jMj ~ Q*Mj m' ' '■^^ fidelity is 

^ JJ' = l/2 Mj = -Jm'j=-J' 



From Eq. 17 we find a lower bound for the fidelity: 

7t K 



Fp,o(E,f/) > 1 



22\gz\ 



(18) 



• If no photon is detected, then measure the photon inside 
the cavity in the rotated basis |x±) and if the outcome 
is Ix-) apply local gate S^^ as in Sec. Ill A. End. 

• If a photon is detected at time t apply the local correc- 
tion gate WcoiT(f )■ Repeat. 

In the above protocol, the unitary U is approximated after a 
sequence of m — 1 clicks and a final null count, an event which 
occurs with probability: 

Pm= dt,„p„a\l{tm) ■ l dt,„^\ . . .dtipi^i{t,„) 
JQ Jo Jo 

Agz\ 2\gz\ 
where the integration measure is dt/t = I /z and 



Pnull(fA;) 



Pdetitk) = 1 -Pnunitk)- 



1. Protocol with a detector 

The above computation of fidelity may be overly pes- 
simistic because one could adopt a strategy where the out- 
put of the cavity is continuously monitored for leakage of a 
photon during the coupling and altering the protocol accord- 
ingly. Consider the situation where we have a perfect pho- 
ton detector outside the cavity that measures the presence 
of a leaked photon. In the case of a null result, the sys- 
tem evolves via a non-Hermitian Hamiltonian: p(f) — ''j^j^^^ 

where i:[p(f)] ~i(yz-ija^a)p{t) + ip{t)(yz + ija^a). The 
operator basis elements then evolve as 



^null 



■(|0)(0| 



\A,J,Mj){A',j',M'j\®-_ 
+e('2«z'^^'^/2)'|0> (1 1 + e(-'2szM/-K/2)f 1 1^ 

_^g(-;2i.z(Mj-M;)-K)r|^^^^|^ 



whereas in the case of detected photon at time t we have 

A^'^^'it) = |A,7,Mj)(A',/,M;|e-'2^z(M.-M;)y^ (19) 

The evolution in the later case is unitary and we can restore 
the system to its prior state by applying the locally generated 
unitary operator Wcorr(f) = e'^^^'^~ . We adopt the following 
protocol to implement U = e^'^^c; 

• Prepare the ancilla as in Sec. Ill A and let the photon 
interact with the system for a time T = K/2gz- 



Assuming for simplicity that the total time to detect a pho- 
ton, apply the correction gate and re-prepare a probe photon 
is roughly T, then the mean fgate and variance Afgate of the time 
to perform the gate is 



^gate — ^ 



.W2|i'z 



eW2|.?z|(7tK/2|^z| + l)-l 



< 2t, 



T;^g7tK/2|^z|(g7tK/2|j;z|(nK/2|^z| - 1) + ^)-^/\gz\ 



g7:K/2|^z|(jtK/2|gz| + l)-l 



< V2t. 



Notice that these values are bounded above. 

The fidelity is easily calculated by noting that the action of 
each map associated with a detector click is a unitary opera- 
tion which is undone by a correction step. So the only step 
in the protocol that acts non trivially is the step with the fi- 
nal null count. Evaluating the fidelity as before but using the 

,Mj.M',, , „ , 

expression A^^^y (f ) we find 



Fp,o(£,f/) = - l+cos2(2x) + - 



2g7CK/4gz 
I gKKl2g2 



sin2(2x) 



2v.2 



= 1- 



Wz 



sin2(2x) + 0((K/2gz)'), 



10 



where the expansion is valid for K/|gz| ^ 1. 

In fact we can do better Using our knowledge of the nec- 
essary gate time T in the event of a null detection it is advan- 
tageous to prepare the initial photon probe state 

(|0)+eW4L?zl|i))/^l+g,:K/2|s.z|, 

Evolution of the basis states is 

^4X^(0 = |A,/,M,)(A',/,M;|0^^-^^i^^(|O)(O| 

^_g(i-2gzM>+7tK/4|gz|-ia/2) | 
_|_g(-i'2«zA/j;+7iK/4|gz|-Kr/2) | 

j^g{-i2gz{Mj-M'j)t+%Kl2\gz\-Ki) 

and A^j'*^''(f) is as in Eq. 19. After a time i a null detection 
gives the target evolution and the fidelity is one. As before, in 
the case of a photodetection event, the system evolution can 
be reversed and the protocol repeated. 

Now the probabilities for a null count or detection of a pho- 
ton are modified to: 

1 _)_gItK/2|gz|-Kft 

P^a\\{tk) = ^^g,iK/2|gz| ' -Pdetlf*:) = 1 - Pnull(fA:) • 

The mean fgate and variance Afgate of the time to perform the 
gate is 

oo 

?gate ~ f^Pm 
m=l 

^ x(l-he"'^/^|gzl)7tK/2|gz| 

7tK/2|gz|+e"'^/2|«zl-l ' 

Afgate — y'f^gate — f^ate 

_ XV(1 +eW2|gzl)(eW2|.?zl(7iK/2gZ - 1) + l)7IK/2|gz| 

~ JtK/2|gz| +e"^/2|A'z| _ 1 

The mean and variance of the gate time is now unbounded 
with increasing K but for K/|gz| ^ 1 the values are com- 
parable to the prior case with the photon probe prepared in 

m+\i))/V2. 

For a situtation with finite detector efficiency r| which can 
be modeled as a rank 2 projector on the photon, the fidelity 
will degrade ultimately to the case of no detector as derived 
above. 



B. Geometric phase gate 

In order to evaluate the effect of cavity decay during the 
the geometric phase gate, we are particularly interested in 
the case where initially A'^-''^j{Q) = \A,J,Mj){A' ,J' ,M'j\ (g) 
\aMj){^M'j\' with |aMj), |Pm)) coherent states. This kind 
of factorization is ti'ue at any stage of spin couphng to the 



field. Using, Eq. 13, the sum becomes an exponential and the 
evolved state is 

p(r) = .^'p(0) 

L e''«aAj,M„A'./',M;|A,/,My>(A',/,M;| 

A,J,Mj,A'J'M'j 

(8> |e-feM+K/2>„^^)(^-feM'+K/2),p^, 1^ (20) 

where 

d{t) = aM.PI,, V,m;(0 - il^Mf + lK'f)^- (21) 

For completeness, an alternate derivation of the evolution us- 
ing the characteristic equation for the joint state is give in Ap- 
pendix B. 

We ignore decay during the displacement stages of the evo- 
lution (i.e. we assume these are done quickly relative to the 
decay rate), and we assume that the system particles do not 
interact with the field during these steps. In order to perform 
logical operations on the protected momory, we do not need 
an ancilla and there are seven time steps beginning with the 
cavity in the vacuum state: 

D(-p)e-"5Vz£)(_oc)e'^3Vz£,(p)g-/x,Vz£)(cc)^ 

as described in Sec. IIIB. Let T5 = T3 = Ti so that the periods 
of spin field coupling are all equal in duration. Notice the 
change in sign of the evolution during the period T3. This 
can be accommodated by changing the sign of the coupling 
parameter gz by e.g. changing the sign of field detuning. In 
order that the field state return to the vacuum at the end of the 
sequence, we choose a'^*^^' , p'^*^"^' . The total sequence then 
yields the output state: 



Pout — 52 '^KJ.Mj.M J' M'j^MjM'j 

KJMjA'J'M'j 

X e'Xsin{.t>+2gzxiMy) |A,7,Mj) (A',/,m;| 

Xe-'Xsin(H2«zXiM;)^|Q^^Q|^ 

where we defined RmjM'j = e''('^'^+'^'^'*^+''^'(') and % = 

lapKe^^""'/^ + /2. This can be interpreted as coher- 

ent evolution with an effective evolution operator 

where T is an effective time for the gate, followed by de- 
phasing in the {A,J,Mj} basis. Matrix elements diagonal in 
are invariant. For m even, the parameters gzii = k/2, 
\a\ = |p|, (|) = ±7i/2, generate 

U - exp[T/(-l)™ |a|25-^(e-3'^"'/2_^e-'"'/2)/2]. 
The strength of the dephasing and decay is then. 
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^m,,m; =exp 



X exp 



{Mj~M'jY + {K/2gzf 

((_l)My _ (_l)M.'/)|a|2e-3W4|gz|(l +eW2L?z|)(K/2g2)2 



{Mj-M'jY + {K/2gzf 



Note that R;^,^ j^j = R*^^ j^, . For K/\gz\ < 1 and Mj ^ M'j, 



(22) 



l-47t|a|2^+47r2(5^)2|a|2(l + 2|a|2) 



Mj — M'j even 



1 -47i:|a|2(l- 



1 



My-M'J 2|gz| 



4|ap((Mj-M;)ii-l)((Mj-M;)(27t|ap+ir)-2|a|-) ^ K ^2 



(5^)2 My odd. 



(23) 



For m odd, the parameters gzXi — k/2, \a\ — |p|, (|) = 0(7t), 
generate 

f/ = expiTiX-l)"^ |a|25-^(e-^^"/2 + e-''"''2)/2]. 

The strength of the dephasing is the same as Eq. 22 but with 
the replacements Mj {2Mj - 1 ) /2 and M'j {2M'j - 1 ) /2. 

As in the case of the single photon mediated gate, the noisy 
implementation of the geometric phase gate commutes with 
the target unitary U — e^^'^^c. The process fidelity is then, 

1 ni/2 J f 

F,,,{'E,U) = -^l^cJc";, ^ ^ 

^ J.J' Mj = -Jm'j=-J' 

where we make use of the fact that /?^/ = R*^^ ^, . The 

magnitude of the coherent state amplitude is chosen to best 
approximate U by fixing |a| = |p| and: 

X = |a|2 (e-3W4|«zl + e-W41szl )/2. 

It is quickly verified that for K = 0, Fpro{'E,U) = 1. By Eq. 23 
we find the lower bound 

F (T. U) > ] 4jr)cK/|gz| 

>'-S(> + S)- 

In closing, note that the process fidelity for the case in 
which the target evolution is unitary can be related to the av- 
erage fidelity via [_o]: 

Fave(2;,f/)= / d\vF{^i\\v)),U\^){\v\U^)) 

^ Fp„{'E,U)D+l 
D+1 

Topologically ordered states have the property that for pure 
states, when the system is divided into two connected do- 
mains, the subsystem entropy scales like the size of the bound- 
ary [ ]. For the surface codes and in the case where the two 
subsystems are just one string of spins and the rest, this im- 
plies that the subsystem entropy of the string is nearly maxi- 
mal because by isotopy the state of any string on the lattice can 
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FIG. 6: Process fidelity of an implementation £ of the many body 
gate U = e^'*^i='"i as a function of cavity decay K (in units of 
the particle field coupling strength gz). Here the time is chosen so 
that g;X = Tz/2 where x is the time spent during each stage of co- 
herent coupling between field and particles. Plots are shown for 
an implementation using the geometric phase gate on m = 9 spins 
(green) and m = 25 spins (orange); and using a single photon for 
m = 9 (red), m = 25 (blue). Also shown are the lower bounds on 
fidelity from Eq. 18 (dashed) and Eq. 24 (dot-dashed). For the 
phase gate, we choose the coherent state amplitudes according to 

is assumed that no decay occurs 
during the field displacement stages. 

be deformed to any other string in the same homology equiva- 
lence class. Hence we expect that the subsystem of spins that 
are acted on during the gate has equal weight on most states 
in its Hilbert space and the measure of fidelity as an average 
measure over pure states is a good one. 



V. FIDELITY WITH OTHER DECOHERENCE 
MECHANISMS 

Up to now we have ignored decoherence mechanisms such 
as radiative decay of the spins into all modes of the electro- 
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magnetic field, and possible sources of noise such as fluctu- 
ating optical trapping fields and stray magnetic fields. Many 
of these effects will be system dependent, however we can 
make some quantitative statements for the case that the noise 
is isotropic. This is a reasonable working assumption because 
in order to obtain the spin lattice models used for protected 
quantum memories, it is assumed that the qubit levels are de- 
generate. Hence absent any special symmetry imposed on the 
environment and control fields, we expect the noise and radia- 
tive decay to act isotropically on the spins. 

A. Collective depolarization 

In the case where the decoherence channels correspond to 
environmental modes that couple coherently to all the spins, 
we can describe the system as undergoing collective decoher- 
ence. This would be case, e.g. for trapped polar molecules 
where the transition microwave wavelength is much larger 
than the optical wavelength spacing between molecules. The 
map describing collective depolarization is described by ap- 
plying a the same random SU{2) rotation to all qubits in the 
system 

£cd(p) = (i-p)P+/' / da[u{£i)f'"p[uH£i)f"' 

jnesu{2) 

J.A.A' ^ 

where p;^^, = (A,y,Mj|p|A',y,Mj). The strength of the 
collective depolarization is parameterized by p such that for 
isotropic decay at rate y, over a time period t, p = 1 — e^^. 
Essentially the collective depolarization erases coherences be- 
tween different J quantum numbers and maximally mixes the 
reduced state within each / block. Accordingly, this map com- 
mutes with the map Eg corresponding to the cavity decay de- 
rived above: E^d ° Ei;(p) = "Eg o Ecd(p), and as before we can 
compute the process fidelity for a process with cavity decay 
and collective depolarization as: 

fproCEcdoEj^f/) {^+\p^,\^+)s,s', (25) 

where "E'lp) = U^'Ecio'Eg{p))U and 

p£' = ^ E \A,J,Mj)s{A',fMj\ 

KJ,Mj.,M J' M'j 
®((l-/')^M,,M;|A,/,My>y(A',/,M;| 

^P^M^ L |A,y,M;)„(A',y,M;|). 

We find 

Fp„CEedO = (1 -p)^pro(E„f/) 

y=(i-(-ir)/4 



The sum can be evaluated in closed form in terms of hyper- 
geometric functions and it quickly decays to zero, e.g. for m 
large the second term scales like )^^^^- • 

B. Independent depolarization 

For environments with a correlation length small compared 
to the lattice spacing, as in the case of optical scattering by 
trapped atoms in an optical lattice we can model the decoher- 
ence as independent isotropic noise on each qubit. Under this 
assumption, the evolution equation will contain an additional 
Liouvillian 

^ k=l a 

where the effective decay rate per particle is Yeff = j 
where y is the spontaneous decay rate and h is the average 
photon number in the cavity. For the single photon mediated 
gate and the geometric phase gate we can assume h ^ \ and 
h— |ap respectively. The corresponding noise map is now 

m J 

■Z^d(p) = (l-mp)p + pETr^[P]®T' 

k=\ ^ 

where Tr^^ is the trace over the A:— th spin and p = 1 — e^^'^f*'. 
Different from the case of collective depolarization, the op- 
erations Eg and Eid do not commute, as can easily be easily 
verified on a two qubit system. 

An estimate for the effect of spontaneous emission can still 
be obtained in the perturbative limit, where mp ~ mjeffX <C 1 
for a gate time x. For both gates this time is T TiA/g^. The 
final state can then be approximated by 

Pout ^ EidoEg(p) 

and the corresponding process fidelity can be lower bounded 

byFp,o(EidoEg,f/) > (l-m79)Fp„(Eg,t/). 

VI. CONCLUSIONS 

We have analyzed the performance of cavity mediated 
many body gates in a spin lattice. To summerize the re- 
quirements for robust gates we require: high quality cavi- 
ties, i.e. low loss rates (K/|gz| ^1), and dispersive coupling 
Y,g <C A. These can be satisfied in the strong coupling limit 
where Y, k ^ 1 . Note that it is possible to have K <C g < Y 
and still satisfy these requirements. Recent experiments [ ] 
reported 3D trapping of Rb atoms in a high finesse opti- 
cal cavity with coupling parameters (g, Y, k)/27I = (16,3,1 .4) 
MHz. Microwave cavities offer the possibility of even bet- 
ter numbers. For example, superconducting strip line cavi- 
ties resonant at microwave transitions frequencies have been 
built ["''] with parameters (g,K)/27i: = (200,0.1) MHz. These 
cavities can be used to trap polar molecules and coherently 
control them on the microwave transitions between rotational 
levels (the linewidths on such excited rotational states are 
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negligible)[^^]. A difficulty here may be interacting with the 
atoms using lasers in the vicinity of the strip line cavities. One 
might try to trap without lasers using self assembly with static 
electric fields but care should be taken to ensure that the under- 
lying lattice model is compatible with the setup. Ultimately, 
we expect the idea of using quantum probes for many body 
control will suggest new strategies for information processing 
in strongly correlated states of matter. 
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APPENDIX A: A SUBSYSTEM CODE 

The Hamiltonian //cp differs from the 2D Ising type 
model introduced in [ ], namely: H' — — /(Lx-iinks^'^^^ + 
Ly- links (^^'^^'''^ + O^'O') + Lz-iinks However, both mod- 
els possess the same subsystem structure. Stabilizer operators 
are generated by the 2(n — 1 ) members of the set {Vj , J 
where the generators are adjacent planes of operators in the 
x — y plane and o~ operators in the y~z plane: 

^ - n ^ = n <^,j^,jmv 

i.k=\ i.i=\ 

The Hamiltonian //cp can encode one qubit of information 

in a subsystem of the total Hilbert space O-ff" ■ For clarity, we 
recall the argument given in [ ] for the subsystem structure of 
the energy eigenspaces. It is understood by considering invari- 
ant subspaces of the Hamiltonian with respect to three sets of 
operators: L, 1^, 'T. The set 'T, which is a group, consists of 
all products of Pauli operators consisting of an even number 
of a* operators in each x — y plane and an even number of (3~ 
operators in each y — z plane. It is generated under multiplica- 
tion by the summands in the Hamiltonian H, i.e. 

(Al) 

The Hamiltonian, in particular, is in the real span of 'T. The 
stabilizer set V, also a group, is generated under multiplica- 
tion as 

1^ = {{Vk .yk}'r=\) (A2) 

V is an abelian subgroup of 'T . Finally, the set L consists of 
operators with an odd number of y — z plane operators Lf — 
n" ^.=1 j ^ and an odd number of i — y plane operators — 

L = {{Lf})/{{Vi^}r{{V^})/{{V?}) (A3) 



This set is clearly not a group (e.g. it has no identity element) 
hvX Lyj'V is. Note that Vf G CTjV e l^^i e L, the following 
commutation relations hold vfv^' — t,t£t^^ — £. 

We can partition the Hilbert space into ±1 eigenspaces of 
the 2(n — 1) independent stabilizer generators {V/, V^^}: 

9{ =®,j(^,.z^^^^ (A4) 

where = (vf , vf , . . . , j ) is an « — 1 bit string of the 
eigenvalues of and = (vf , vf , . . . , j ) is an n — 1 bit 
string of the eigenvalues of V^. Because the Hamiltonian 
^^cp € spangT, its eigenspaces are block diagonal in {v^,y^}. 
Furthermore, because all elements of L commute with ele- 
ments of T, we further decompose the eigenspaces as 

i?<;^,„z = i?<5,z®i^,x,,z (A5) 

Operators in L commute with T and 1^ so they leave those 
spaces invariant. In a given stabilizer eigenspace, any oper- 
ator in L can be reduced to the simple product of one plane 
operator and one plane operator Lf. Because n is odd these 
operators anticommute: {Lf ,Lf } = 0, hence they form a rep- 
resentation of a two dimensional Clifford algebra. By dimen- 
sion counting then: dimi?/^ = 2" ^jj^j (jimj?/^ — 2. 
It is in the subspace L that a logical qubit can be stored. Fur- 
thermore, the logical operators on the qubit subspace corre- 
spond to single plane operators. Writing if if = iLf, the al- 
gebra spanj{{Lf ,Lf ,Lf } then forms a representation of the 
algebra su(2); i.e. they are the logical qubit operators. 

There is another way to see the action of these op- 
erators on the ground states of //cp. Note that //cp is 
time reversal symmetric and the number of spin— 1/2 par- 
ticles in the system is rr' (odd). Hence by Kramer's The- 
orem, each eigenspace has degeneracy which is a multi- 
ple of 2. The eigenstates come in pairs (|A,),LS|A-)) where 
the anti-linear time reversal operator acts as 13\X) — KCl'k.) 
where C is the complex conjugation operation and K = 
11" j k=i{^^'^] j k)- Because the Hamiltonian is real, the 
eigenstates can be chosen real such that any pair are given 
by {\X),K\X)). Now K = --Wj,k=i<^l,MlY;,.k=i<^ij.k = 
-[nti ^^^Vf,]Lf[Y^l-'^^^Vi]Ll But in a given stabi- 
lizer eigenspace !J-([,x ,,z, the action of this operation is K = 

^iIl["kJi^^ ^i^2k)^^^^ ■ Restricting to the subspace vf = 
vf = IV/, fe, we have K — —iL\. This then defines (— ;) times 
the logical Y operation on that subspace. From the commuta- 
tion relations the operators Lf and if are the logical X and Z 
operations respectively. 

A logical CNOT operation can be done transversally be- 
tween two code blocks. This follows by considering the action 
of the CNOT, j operation on the Pauli operators: 

CNOT,;yX,CNOTi. , = XiXj , CNOT,;yZ/CNOT, j- = Z, , 
CNOT, jXjCNOTi j = Xj , CNOT, ^ZyCNOT, = Z,Z,- , 

(A6) 

Hence, by the group homomorphism, the joint stabilizer 
groups for the control / and target J logical qubits V y. V is 
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preserved under conjugation by CNOT®" = CNOT/ j. This 
is easily checked by noting that even numbers of planar , 
operators get mapped to even numbers of planar operators of 
the same type. However the logical operators are acted upon 
nontrivially. Specifically, we have the same relations as in Eq. 
A6 but with logical operators replacing the physical qubit op- 
erators. Hence the transversal CNOT is a logical CNOT on 
the code blocks. Given the ability to generate arbitrary rota- 
tions e'*^ , e'*^ and the CNOT gate, and measurements of L^, 
exactly universal quantum computation is allowed. 

The corresponding code is a [[rr", l,n]] code, i.e. it encodes 
1 logical qubit in rr' physical qubits with a distance d = n. 
This code can detect up to — 1 arbitrary errors and correct 
\{n— l)/2] errors which is the maximal length of a an ar- 
bitrary error string with an unambiguous action. Error cor- 
rection is done by finding the minimum Hamming weight n 
bit string consistent with the stabilizer measurements and ap- 
plying single spin(phase) flips on those planes correspond- 
ing error locations. An example of a worst case er- 
ror which saturates these numbers is the length £ error string 
E — Y]Li'=il '^i j k- error string has two boundaries which 
flips the sign of the stabilizers measurements and V^^_|_f- 
For t < \{n— l)/2], the error string creates two boundaries 
which are detected and appropriately corrected. But for longer 
strings the error correction proceedure implements a logical 
error on the code. 

Some requirements for robust information processing in the 
above code are: 

• Preparation of the system in the ground subspace of H. 
Presumably can be accomplished by cooling the sys- 
tem to a pure separable state of a local Hamiltonian 
Hq = YJi j k=\'^] j k ^^^^ adiabatically turning on H. It 
would be necessary to check that the adiabaticity re- 
quirement was satisfied by estimating the gap of the 
time dependent Hamiltonian H{s) = {\ —s)Hq + sH,s = 

t/Te [0,1]. 

• Projection onto a fiducial logical state in the ground 
subspace. This can be done by measuring mz = (if) 
and assigning the logical state |0l)(|1l)) to outcome 
mz = ±1. Such a process also allows measurement of 
the state in the logical Z basis. 

• Encoding a quantum state and performing single qubit 
rotations. In order to do so it is necessary to be able 
to implement single logical qubit operations. To gener- 
ate the continuous group SU(2), the continuous gate set 
{e'^l ,e'^^} suffices. One may rather demand only a 
discrete gate set that generates a group dense in SU(2). 
One such gate library that can be done fault tolerantly 
is {L'' = e'l^fe'l^re-'l^i ,CNOT}. Since the is 
a product of & operators on all qubits, single logical 
qubit gates could be performing by emersing the crystal 
in one direction into the cavity, first one plane to gener- 
ate gates via Lf and then the entire crystal to generated 
gates from L^. The CNOT can be performed transver- 
sally between two code blocks by performing physical 



CNOT gates in parallel with control physical qubits 
in one code block acting on rr' target physical qubits in 
the target code block. This one step parallel operation is 
difficult to do with local operations, however one could 
perform the gate locally by performing parallel CNOT 
gates between the physical qubits in the bottom most 
x~y plane of the control logical qubit and the target 
physical qubits in the top most x — y plane of the target 
logical qubit. A series of « — 1 such steps where the 
X — z planes of the logical qubits are cyclically shifted 
realizes the logical CNOT. Each cyclic shifts can be 
done in a linear number of parallel planar SWAP gates. 

• Measuring stabilizer operators is necessary to detect 
and ultimately correct errors. This demands measur- 
ing the set of 2(n — 1) independent stabilizer generators 
{V,^ , V^} which are nearest neighbor planes of products 
of all or all a' operations. 



APPENDIX B: ALTERNATE DERIVATION OF DYNAMICS 
DURING THE GEOMETRIC PHASE GATE 

The evolution of the joint system of spins and field in a co- 
herent state basis was derived in Sec. IV B by integrating the 
equation of motion including cavity decay. Here we provide 
an alternative derivation using characteristic equation for the 
state. We begin by transforming the evolution in Eq. 12 to an 
interaction picture via 

p/(f) = e'^^'p{t)e-'^^', a,{t) = e'^^'a{t)e-'^^' , 

we have 

P/(0 = K(fl/P/flj' - ^a]aipi - ^p/flj"fl/). 

Now fl/(f) — e^'^^^''~'a hence, 

Piit) = K(ae-'2sz^''p/e'2«^-'''a+ - ^a^api - ^pia^a). 

In the interaction picture: 

Af (t) - |A,7,M,)(A',/,M;|(f)® |aM,(f))(pM;(f)l, 

where: = aM,(f)e'2.?z'^', p^, (f) = p^/ (Oe'^i'^'^^), To 

derive the evolution during decay we use the characteristic 
function 

X{t)^Tr,[Af''^'^{t)e'^'e-^'"], 
(where the trace is taken over the field) such that 

X(f)=TrHAf"''^(0e'"'e-'*1 
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Using the relations: 
we obtain 



dX*dX 2 \ dX* 
(^^M,{t)X-aMjit)X*y. (Bl) 



K /s* 

2' 



The equation of motion for the operator X can be solved by 
the method of characteristics. For the states of interest here 
we make the ansatz: 

X{t) ^Cit)e-^'''^"j'''^e^^"j^'\A,J,Mj){A' ,J' ,M'j\. (B2) 

Evaluating the time derivative and setting this equal to Eq. B 1 
we find 



(P^,e-'°/2|aM,e-^/2)|A,7,M,)(A',/,M;|. 



Hence, 



Jo 

+ Ke-«(l-e-'2«^'(^^-0)]}, 



where we have chosen the integration constant so that 
c(0) = 0. Notice that for gz — 0, i.e. pure decay, then 
c{t) — 0. This should be the case as then |aMj)(pM'l ^ 



(PM^IaM,)'^-'" )|aM,*?-^/2)(p^,g-K,/2|^ Furthermore, for 

Mj = M'j then c(f ) ~ 0. To account for all phases and decay 
we introduce: 



d{t) = c(f) 



ccm, 



Finally we arrive at 



p(0 = .^'p(o) 

KJMjA'J'M'j 



which is the same evolution as derived in Eq. IV B. 
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